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Abstract 

We discuss the universal behaviour of scattering cross sections in the hmit of infinite 
rapidity separation between all produced particles, and illustrate the behaviour explicitly 
for the production of n jets, W + n jets, Z + n jets for n = 2,3, 4, and for H + 2,3 jets. 
We give a set of rules for constructing scattering cross sections, which are exact in the 
given limit, and order-by-order reproduce well the full fixed order results when applied to 
LHC phenomenology. The approximation includes both real and virtual corrections, and 
is sufficiently simple to allow for the regulated all-order perturbative sum to be explicitly 
constructed. We test the expected accuracy by comparing results ordcr-by-order with full 
fixed-order perturbation theory for the processes mentioned above. 
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1 Introduction 



Achieving the fuh potential of the LHC in furthering our understanding of particle physics will 
challenge our understanding and description of events with multiple jets. Such events arise 
both within the Standard Model (SM), and in many models of physics beyond the SM. Apart 
from furthering our understanding of the solutions to the SM equations of motion, a better 
understanding and description of the multi-jet predictions arising from the SM is necessary in 
order to fully disentangle this contribution from that which might arise from outside the SM. 

The true complication of an observed jet in terms of its constituent hadrons can currently 
only be described within the context of a "General Purpose Monte Carlo", implementing a 
parton shower and hadronisation model as described in e.g. pbiSJ. The parton shower is based 
on a resummation of soft and collinear radiation. However, while obtaining a good description 
of the structure of each jet, this description severely underestimates the rate and hardness {p±- 
spectrum) of multi-jet samples. Better predictions of such quantities in exclusive (meaning 
fixed number) jet samples have so far been obtained in fixed order calculations; for processes 
allowing more than two light partons in the final state, such predictions are currently limited 
to NLO accuracy. However, putting aside issues of hadronisation, many other questions which 
are important for the LHC programme still cannot be answered satisfactorily in the relatively 
simple description of the final states obtained at such low orders in the perturbative expansion. 
Such questions range from the effects of central jet vetos in samples of Higgs boson plus jets, to 
the simple question of the transverse momentum spectrum of weak gauge bosons (both at large 
and small transverse momenta) or the relative weight of various jet multiplicities in inclusive 
samples ("inclusive" here meaning sum over any number of jets). 

A combination of fixed order calculations and further parton shower resummation has been 
achieved, which so far corrects the approximation obtained in the parton shower to full tree-level 
accuracy (for a fixed jet multiplicity IMM) or NLO accuracy [HIT] (so far in processes which at 
LO have just zero or one light parton in the final state). Such matching schemes ensure that 
a given process is described correctly at least at (N)LO, while all higher order corrections are 
estimated from the parton shower. 

All resummation schemes are built upon a particular kinematic limit in which the pertur- 
bative corrections simplify, allowing for all-order approximations to the full perturbative series. 
For the parton shower, this is the soft and the collinear limit, i.e. emissions under small invariant 
mass, and the building blocks are e.g. parton splitting functions. In the current study, we will 
be guided by the simplifications of the perturbative corrections in the case of large invariant 
mass between all emissions. We will seek to catch the part of the perturbative corrections which 
controls additional jet production, but ignore the collinear behaviour which gives rise to the jet- 
sub-structure. We will obtain an approximation to both virtual and real-emission corrections, 
which allows us to build an approximation to the regularised, all-orders matrix element for each 
exclusive (resolved) parton multiplicity. The focus of this paper is to obtain the relevant uni- 
versal building blocks for multi-jet predictions. Furthermore, we will investigate the accuracy of 
this approach by comparing the lowest order predictions obtained in our simplified approach to 
the full tree-level QCD perturbative series for n— jet-production, (n = 2,3,4), both "pure" and 
in association with a W,Z or Higgs-boson. 

The approximate or "f-channel factorised" matrix element will be sufficiently fast to evaluate 
numerically that the resummation can be constructed by explicit summation over the exclusive 
final states, and thus any event analysis can be performed by simply imposing jet algorithms 
etc. on the exclusive final states ("exclusive" meaning fully differential in the momenta of all 
produced particles). We will leave the phenomenological implications of the resummation of 
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each of the four processes discussed here to future studies. 

In Section [2] we wih introduce the Multi- Reg ge- Kinematic-limit, which hes at the heart of 
our approach, and discuss and illustrate the universal behaviour of scattering amplitudes in 
this limit. In Section [3] we will construct all the building blocks for the relevant amplitudes. 
In Section m we will check process-by-process, multiplicity-by-multiplicity the level of accuracy 
at which the approach reproduces the fixed order perturbative expansion in terms of not just 
cross sections but also differential distributions, where these can be obtained with standard tools 
implementing the full tree-level processes. By doing so, we hope to instill trust in the all-order 
approximation which can be built using the elements described in this paper. 



2 The High Energy Limit of Scattering Amphtudes 

We will start by studying the scattering matrix elements of a gluonic 2 — )■ n, n > 2, process 
in the so-called Multi-Reg ge- Kinematic (MRK) limit, where the scattering momenta pa,Pb 
pi, . . . ,Pn in terms of transverse momenta and rapidity y = In ^ E-p'' ) ^^^^^ the following con- 
ditions 

Vz E {2,. . . ,n - 1} : > yi > j/i+i 
Vi,i : \pi±\ ~ \pj±\, 

or alternatively 

Vi,j : \Pi±\ ^ \Pj±\,Sij oo, (2) 

where Sij = 1 Pi-Pj and s = 2 pa-Pb- The notation yi ^ yj really means yi — yj — S" oo. The 
transverse directions are with respect to the incoming partons pa,Pb (i-e. transverse to the beam 
line), and the limit requires the transverse components to be kept fixed (i.e. not growing with 
s) as Sij, \yi — yj\ ^ oo. We can take pA to be the direction of positive light-cone momentum. 

Explicit calculations of up to 2 —t- 4 gluon scattering [8], and a dispersive analysis of the 
2 — )■ n process [9,10j showed that in the MRK limit, the 2 — >• n scattering amplitude is dominated 
by the behaviour dictated by the poles from t-channel gluon exchanges, t-channel here refers to 
the picture which arises of a string of effective vertices for on-shell gluon production, which are 
connected by off-shell (t-channel) gluon currents, see Fig. [TJ The t-channel momenta are here 
defined as 

i 

ti=PA-^Pj- (3) 
i=i 

In Section [3] we will illustrate how the following result arises, and re-derive the form of the 
building blocks of the approximation. Here we will just quote the result for the LO colour and 
spin summed and averaged scattering matrix element for 2 — t- n on-shell gluons in the MRK 
limit: 



M 



[j-J; \Pi±\^ ) \Pni. 



One notices that obviously, after the MRK limit is taken, there is no dependence left on the 
rapidities of the gluons (or rather, the dependence on the rapidities is lost). The dependence on 
the centre-of-mass energy s in Eq. Q is left in, although of course in the MRK limit, s — )• oo. 
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Figure 1: Picture of effective vertices connected by t-channel exchanges 



The partonic cross section is found by dividing the square of the scattering matrix element by 
the flux factor, which is proportional to s^, so the MRK limit of the partonic cross section is a 
very simple function of only the transverse momenta of the emitted gluons. 

The results of Ref. [OHllj extend beyond just pure gluon scattering; the picture which arises is 
one in which the description of a given 2 — )• n scattering process in the MRK limit factorises into 
a product of effective vertices for particle emission, with the vertices connected by propagators 
according to the ordering in rapidity of the emitted particles. One immediate result is that in 
the MRK limit, the 2 — )• n scattering process is dominated by rapidity orderings which allow 
pure gluon exchanges (or more precisely: exchanges of the particle of highest spin) between the 
scattering vertices. 

Another immediate result of the factorisation (arising from the assumption of a hierarchy in 
the components of plus- and minus-momenta) is that in the MRK limit, the scattering matrix 
element for e.g. gg — )■ ggg, qg — )■ qgg and qQ — )• qgQ differ only by colour factors - this is a 
generalisation of the results leading to the "effective PDF" -approach to 2 — )• 2 scattering |12j . 
The results equivalent to Eq. Q for the qQ and g^^-initiated scattering processes in the MRK 
limit are then 



-r-^MRK 



-t-tMRK 



g^CF 



NS-i 

4s^ g^ Cf 





g^CA 


ImP ) 


\Pn±\^ 


^g^CA) 


g^ Cf 




\Pn±\^ 



(5) 
(6) 



with the additional information that all final states not ordered in rapidity according to the 
indication in the subscripts are suppressed (by powers of Sij). The difference in colour factor 
when gluon 1 or n is replaced by a quark is Cf/Ca, since the difference in the summed and 
averaged colour factor is the replacement of ^''^"^''/(^^-l) = CAd'"^/{N^-l) with %tjjNc = 
5'^'^ /{2Nc). Eqs. arise also as the expansions to fixed order of the solution to the BFKL 

equation to leading logarithmic accuracy. 

The results of Eqs. Q-dU]) are summarised in Fig. [51 which illustrates the behaviour of the 
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Figure 2: The scattering matrix elements for (a) qQ — )• qQ, qg — )• qg and gg — )• gg, (b) qQ — )• qgQ, 
qg qgg and gg ggg and (c) qQ qggQ, qQ qggQ and gg gggg, as a function of 
the rapidity difference A between the partons, defined in Eq. ([S]). The MRK limit is the value 
obtained from Eq. ([6]). 
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The full tree- level matrix elements are extracted from MadGraph [13] , and the momenta of the 
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final state particles are chosen as 

Pi = {k± cosh(yj), k± cos(0i), k± sm{(j)i), k± sinh(?/j)), k± = 40GeV, 

(a) : yi = A, 7/2 = -A, 

(b) : yi = A, 7/2 = 0, ys = 

(c) : yi = A, y2 = —, ys 

The same universal behaviour is seen in processes where a W, Z 01 H boson is produced in 
association with jets, as shown in Figs. [SHSl The momentum configurations for these plots are 
given in Appendix |Bj The rapidities of the jets in each case are as in Eq. ([8]), except for the 
H + 3j process, where the jet rapidities are A, —A/3 and —A respectively. The MRK limit for 
the W processes is taken from Ref . [H] while that for the Higgs processes is taken from Ref . |15) . 

As can be seen on the plots, the correct MRK limit is obtained for all the channels. However, 
while the simple formulae of Eqs. (Il])-(l6]) do indeed describe the MRK limit of the full amplitudes, 
this limit is approached only outside the region of relevance for the LHC. The power vested in 
this implementation of the factorised picture in terms of ability to calculate a specific limit of 
the n-gluon amplitude is turning out not to be relevant for LHC phenomenology, if implemented 
according to Eq. Q or the BFKL equation. 

In contrast, the universality of the qQ, qg and gg-channels, which arises naturally within and 
to some extent even implies the picture indicated in Fig. [U is reached in the region relevant for 
the LHC phenomenology. In fact, the qQ and qg channels behave very similarly over a very wide 
range of rapidities. This poses the question whether one can construct better building blocks 
still within a factorised picture, which would allow one to get estimates for the all-order cross 
sections, but crucially in a wider kinematic region which would be relevant for the LHC. This 
is the focus of the next section. 

3 The Method Of Currents 

The aim of this section is to construct a method for approximating the hard scattering matrix 
element for jet production (pure or in association with a W/Z/ H-hoson) to any order in a^. 
The method should have the following characteristics :- 

1. Inclusivenes: The approach should capture both real and virtual corrections to all orders 
for a given process. 

2. Simplicity: 

(a) All orders exclusively: In order to allow for arbitrary analyses of the process, the 
all-order result must be constructed as an explicit sum over n-particle final states, 
with access to the momenta of all emitted particles. This obviously requires that 
the evaluation of the scattering amplitude for any number (necessary) of particles 
is sufficiently fast to allow for the n-body phase space integration to be performed 
explicitly. 

(b) Cancellation of IR poles: The formalism has to be sufficiently simple that the cancel- 
lation of IR poles between real and virtual corrections can be organised while keeping 
the all-order summation simple. 

3. Accuracy: 
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Figure 3: The scattering matrix elements for jet production in association with a W boson, for 
(a) 2, (b) 3 and (c) 4 jets. 




Figure 4: The scattering matrix elements for jet production in association with a Z boson, for 
(a) 2, (b) 3 and (c) 4 jets. 




Figure 5: The scattering matrix elements for jet production in association with a H boson, for 
(a) 2 and (b) 3 jets. 
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The obtained results must reproduce the full perturbative result order by order in 
in the limit of infinite invariant mass between all partons. 

At the same time, the result must maintain relevance in the kinematic regime of 
TeV-scale colliders; we will calculate the cross section and kinematic distributions 
for 2, 3, 4-jet production alone and in association with a W, Z or H boson within the 
simple procedure allowing an all-order construction, and compare the results to those 
obtained with full, fixed-order calculations. 

Our approach will be based on capturing the behaviour dictated by poles in the t-channel 
momenta of the scattering amplitude (as defined in Eq. ([3]), where the momenta pi are rapid- 
ity ordered). Furthermore, in this paper we will focus on the rapidity ordering of particles 
which allows colour octet exchanges between all neighbouring (in rapidity) particles. It is well- 
known |lU y i6 t [T7] that for a given 2 — )• n process, these rapidity-orderings will form the leading 
contribution to n-jet production, and all other orderings will be suppressed by powers of the 
invariant mass between jets of flavours which cannot be connected by a colour octet exchange. 
The relevance of this property of the scattering amplitudes for inclusive cross sections (i.e. not 
just as an asymptotic argument) was investigated for Higgs boson production in association with 
at least two jets in Ref. |18pi9j. 



3.1 Current-Current Scattering 

Based on the close resemblance between the various partonic channels in2— )-2,2— )-3,2— )-4 
scattering, we choose to model all channels on qQ — )• q(g . . . (5r)Q-scattering. The similarity in the 
behaviour of the various partonic channels in the MRK limit, where the t-channel dominance 
has set in, is further supported by the behaviour of the three-gluon vertex and the quark current 
(Fig. [6]) in the limit ~ pi (using the standard spinor formalism): 

{IHA) ^ 2^:^/1 

BA^^el^AiPA+Pir 9"^^' + {q-Pir^ 9"''' + {-q-PAY^Q^n ^ ea ■ el 2p\,^. 

The building blocks of the traditional (B)FKL formalism are termed impact factors (describ- 
ing the coupling of e.g. two on-shell partons to an off-shell (t-channel) gluon) and the kernel 
(describing the evolution of the off-shell gluon under the emission of partons). These building 
blocks can be obtained by taking the MRK limit of specific processes, see e.g. Ref. [2D] and refer- 
ences therein. However, with building blocks obtained in this way, the resulting approximation 
for the scattering amplitude will have lost all dependence on the rapidity difference between 
particles (since this has been taken to infinity), and only the limits in Fig. [2] are obtained. When 
the BFKL equation is used, this approximation is then applied in all of phase space, resulting 
in a gross over-estimate order-by-order of the scattering amplitude. 

Instead of relying on a kinematic limit to ensure the t-channel dominance of the scattering, 
we will instead obtain the building blocks from the process qQ — t- qQ, which consists only of a 



PA,l^A 




PA 



Pi 



V u 
Figure 6: The three-gluon vertex and the quark current shown are identical in the limit pA 
Eq. ©. 
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t-channel gluon exchange. This means that the lowest order results for the qQ — )• qQ process will 
be identical in full QCD and with the effective rules. The benefits are even greater for e.g. W or 
Z plus jets, where in the standard approach (see e.g. Ref. |T3])) kinematical assumptions on also 
the weak boson (and its decay products) as well as the jets had to be applied in order to obtain 
the building blocks for a resummation. We will see that the results for Higgs boson production 
in association with jets |15 y i8 1 [T9] can also be improved using this new formalism. 

When the effective Feynman rules we derive here are used in constructing all-order results, 
the differences between the various partonic channels at small rapidities (and the corrections to 
the approximations for multi-parton production) will be taken into account order-by-order by 
matching to the full fixed-order perturbative result using the same approach as in Ref. |18yi9j. 

In the spinor notation, the (colour and coupling stripped) matrix element for the process 
qQ — )■ qQ for negative helicity quarks reads 



(l|^|a)Y(2|^l^) 
2 [a b] (2 1) 



(10) 
(11) 



We will use the scattering of quark currents as the basis for our framework, since these consist 
of t-channel gluon exchanges only, and explicitly exhibit the factorisation into two components 
(spinor strings), each depending only on the momenta along each quark line. This is also 
obviously true for all other helicity configurations. 

Let us denote the spinor string (for helicities ha,hi,hb,h2 of the quarks) appearing in the 
amplitude as 



(1 hMa ha) g^" (2 h2Hbhb). 



(12) 



This complex number will be calculated using an explicit representation for the spinors (see 
Appendix |A|) . and we will denote the sum over helicities of the absolute square of this number 
by 



\S. 



qQ—^qQ\ 



s: 



haht^hih'z 
qQ^qQ 



(13) 



Of course in this case non-zero contributions arise only when ha = hi and hb = h2- 

The colour and helicity summed and averaged matrix element for the scattering process 
qQ — )• qQ is then 



M 



qQ^qQ 



1 



4 (iV2 



c 



1) 



\S, 



qQ-^qQ\ 



9^ Cf f 
5' Cf 



(14) 



with ti = {pa — Pi)^ and t2 = {—Pb + ^2)^ (^i = ^2 in this case of a 2 — 2-process). The 
superscript t is meant to indicate that this is the "t-channel factorised" -approximation, but this 
is of course exact in the channel qQ — )• qQ. The approximations for the gluon channels are 
obtained by multiplying Ca/Cf for each pair of quarks replaced by gluons. Fig. [2] shows that 
this is an extremely good approximation. 
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The spinor formalism clearly displays the factorisation in the t-channel of the scattering, 
i.e. within the spinor formalism, the "impact factors" are clearly identified as the quark currents; 
in contrast, the "spinor product" rewriting of Eq. (jlip mixes momenta from the two quark lines. 
The standard procedure for extracting impact factors using the helicity formalism [5D] applies 
the kinematic approximations valid in the MRK limit. In terms of invariants, the square of the 
colour and spin averaged and summed scattering matrix element for qQ — )■ qQ is 

9 9 (15) 

The li^-terms arise from scattering of quark currents of different helicities, and spoil the factori- 
sation implied in Eq. ([6]) from being exact, necessitating the consideration of kinematic limits of 
the squared scattering matrix element. This despite the fact that for processes which proceed 
only through a t-channel gluon exchange, the starting expression in terms of spinor strings is 
already factorised. 

In the MRK limit (of infinite rapidity separation between the scattered partons) , the expres- 
sion for the colour and helicity summed and averaged matrix element simplifies to the 2-jet part 
of Eq. ([6]) since all allowed helicity scatterings give the same result in the limit and 

|[a6](21)| = s, t^-\p±\\ (16) 

This is the lowest order results in Eqs. dH)-®. 

3.2 Multi-Parton Production 

In this section we will develop the picture of the scattering of two quark currents to take into 
account the emission of additional gluons. We first consider adding one extra gluon to the 
qQ — )■ qQ scattering we have taken as our model so far; this may be emitted from the t-channel 
gluon or from each of the external quark lines. Fig. [71 

Firstly the three gluon vertex in the t-channel emission gives a contribution of 

tlt2 

(17) 



A = -rr^ uipih^PA) u{p,Mpb) e; ((gi + q2rg^'' + {P2 - 92)^9"' - (.qi+P2r9^n 



where Cg = T^^^^^f^^^'^T^.^^^. In the MRK limit, we can use Eq. ([9]) for the spinor strings, and 
qi = P3 - PB + P2 and q2=Pa-Pi- P2 to get 

Ag ^ e; (-p^(s3B + 2s2b) + P^(siA + 2s2a) + (gi + q2rs) . (18) 

The MRK limit (Eq. ([T])) gives S2b ^ S3_B and S2a ^ siA so we are left with 

^ 2p'/-f + 2p^/-f + iq, + q2r). (19) 
PA Pi 

PB P3 PB P3 

Figure 7: We add contributions where the 3rd jet (red) is emitted from the i-channel gluon, and 
from each of the four external quark lines. 
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We now add the contributions from the emissions of gluons from each external quark hue. We 
treat these emissions as soft, and use Eikonal factors which are vahd for a gluon emitted between 
the quark jets in the MRK hmit. These add to give 



Aq = AqQ^qQ X (igs) (^C 

The d are the relevant colour factors: 
Now in the MRK limit, pA 

SqQ~^qQ 



p 

Pi 

Pi ■P2 
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PA ■P2 



+ C3 
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P3 ■P2 
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Pl,PB-^ P3 and AqQ^qQ = SgQ^qq/t, SO 



An 



{-igt)el \{Ci-CA)qi 



p 
Pa 



PA ■P2 



+ (C3 - Cb) ql 



p'b 



PB ■P2 



(20) 



(21) 



(22) 



However, (Ci — Ca) = iCg and (C3 — Cb) = —^Cg so that all five contributions give the same 
overall colour factor. (The same pairing up is seen when a quark line is replaced by a gluon line, 
only here it is the Jacobi Identity which can be used to give a single overall factor.) 

We could then add Eq. ([22]) to Eq. (fT9]) : this gives the form used in [18l[19]. However we 
choose to first reinstate the full kinematic structure of Eq. ()20p in order to capture more of the 
original process. We also adapt Eq. (fT9|) to take account of both {j)a-,Pb^ and {^1,^3} to give 
the following approximation for the qQ — t- (/(/Q-scattering matrix element 



AqQ^qqQ = 9s Cg £* ''^^ J^^'^ ^^(91,52) 



(23) 



where 



y^{q\.q2) = -{qi^q2) 



p 

,Pa 

2 

Pb 



_^ P2-PB ^P2-P3 



P2 - PA PA- PB 



4 



^ P2 - PA _^ 



PA-P3 
P2 -Pl 



P2 -PB PB- PA PB - Pl 



+ PA^ Pl 



PB ^ P3- 



(24) 



This form of the effective vertex is gauge invariant; the Ward Identity, pg ■ V = can easily be 
checked. The spinor structure is exactly as for qQ — )• qQ and so we find 



M 



qQ^qgQ 



1 



4(Ag 



1) 



\\SqQ^qQ\ 



(25) 



-9'Ca 

tlt2 



V^{qi,q2)V^{qi.q2. 



The results of this formalism are compared to the full tree-level matrix element for qQ 
in Fig. [H^a). The results from the two methods are indistinguishable by eye. 
It can be checked that in the MRK limit, Eq. (p^ gives 



qgQ 



-V>'{qi,q2)V^{quq2) 



'2± 



(26) 
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Figure 8: Results for the matrix elements of (a) qQ — )• qgQ and (b) qQ — )• qggQ processes 
respectively. The momentum configurations are as for Fig. EJ^b) and (c) respectively. 



The l^jxp terms in the numerator cancel with the factors of tj in Eq. (|25p to give an MRK limit 
which agrees with Eq. ([6]). Motivated by the structure of Eq. ([6|), we describe a general process 
with n jets in the final state (Fig. [1]) as 



M 



qQ^qg...gQ 



1 



4{m 



-yy \\SqQ^qQ\ 



9^ Cf 



9^ Cf 



n-2 

n 

1=1 



-9'Ca 



V^{qi,qi+i)V^{qi,qi+i, 



where V^{qi,qi-^-i) is the obvious generalisation of Eq. ([2 
VP{qi,qi+i)-- 



{qi + qt+iY 

,2 



p 

Pb 



Pi+l ■ PA 



_^ Pi+l ■ PB _^ Pi+l ■ Pr 



PA-PB 
^ Pi+l ■ PA ^ 



PA -Pn 
Pi+l ■ Pi 



Pi+l -pB PB- PA PB ■ Pi 



+ PA'^ Pi 
-pB ^ Pn- 



(27) 



(28) 



Each emission vertex comes with a factor of in |A^p because (by the discussion above) there 
is a j'j-(j'^''cjabc^ place of a 6^'^ 6^'^ compared to the process with one fewer jet. 

The results for 4 jet final states are shown in Fig. El^b) compared to the full matrix element, 
and show that the formalism performs well. A more detailed study of the results obtained in 
this formalism for kinematic distributions is presented in Section [4.11 



3.3 W and Z Boson Production in Association with Jets 

To apply this formalism to W and Z production with jets, we follow the same structure as the 
pure jets case: modelling both qQ and qg channels on qQ -+ q{g . . . g)Q + {V -+)ii. When a W 
or Z is emitted from a quark line, the two diagrams shown in Fig. [U] need to be considered. The 
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Figure 9: The two diagrams which contribute when a W or Z boson is emitted from a quark 
current Une joined to the rest of the diagram. 



gkion hne with free index // hnks this current to the rest of the diagram. The two diagrams add 
to give 



1 ; I Ui'JaV£- 



t 



oU 



t 



iU 



(29) 



We can rewrite this in terms of manageable spinor-matrix-spinor pieces as: 



J^ii,i,l,o) 



+ 



2pf[uo-ff'PoUi\ - [uojt'PoUe][ua''Pm] - [uoj''PoUg][ugj''PiUi]\ (30) 



X [u^aPeu^] 



where Px is the projection operator P± = (1 ±75)72 according to the hehcity of particle x. The 
current is zero unless Pq = Pi- 

There are two classes of subprocess to consider, those where 

1. the W or Z boson can only be emitted from one quark line, which is the case in all qg 
channels and certain qQ channels for W production, or 

2. the W or Z boson can be emitted from either quark line. 

In the following, without loss of generality, we pick a specific example subprocess for each case: 
for case 1, ud — )• dd{W — )-)f^^ and for case 2, ud — )• ud{Z -^)M. 
Case 1 is straightforward and the analogue of equation (|12p is 



^-/l(,->-/l2- 
ud—^dduii 



Jv^{ijJ,o) (2 h2W\bhb). 



(31) 



Again we use double vertical line notation to symbolise the sum over helicities of each jS'p. We 
write 



\S, 



ud^dup 



|2 _ 9w 



M 



ud—^duiM 



-ht^-h2 12 



— \S' 



(32) 



1 



4 (Ag - 1 



\s., 



ud^duild I 



Cf 



(33) 
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where now ti = {pa — Pi — Pu — Pe)'^ and t2 = {—Pb + P2Y (still ti = t2 in this 2j case). 

In case 2, it is similar but slightly more complicated as one must divide by ts before squaring 
so we get 



M 



1 



4 (Ag - 1) 



{g^ Cf) ■ [g^ Cf) 



E 



ha, hb, hi, 
h2, hi, hi 



ku Jvf.ia,e,ll) 5^^" (2 h2Hbh)-- 



+kd (1 hi\v\a ha) g^'^Jvuib, 1,1,2) ■ - 

tb 



(34) 



where ku^d are the different couplings of top and bottom emission, k^u is the Z coupling to the 
leptons, tt = {pa -Pi-P£- Pif and tb = (pb - P2 - Pi - Plf' ■ 

As in the case of pure jets, these results are exact for qQ — )• qQ + U. We now want to extend 
this to 2j — )• nj + U for n > 2. In case 1, the inclusion of further gluon emissions (ordered in 
rapidity) is done in exactly the same way as for the pure jet case, outlined in section 13. 2[ The 
results are shown in Fig. [10] for W + 3j and W + Aj. 

Case 2 is very similar, except the Lipatov vertices must be included in the analogue of 
equation ()34p before squaring. For example, consider our sample ud process with n jets in the 
final state. Fig. [TTJ To make the (simple) structure clear, we introduce the shorthands 



at = Jvf.{a, 1,1,1) gf"" (2 h2Hb h) 



and 



afe = {lhi\u\a ha) g'"'Jvu{b,i,i,2) (35) 



for the current parts. Then the | • • • | expression in equation (j34p becomes 



ku at 



Vf>^iqt„qt,)...VP-^iqt„_„qt^_,) , , _ V^^ {qb„ qb,) . . .V^" 



tti--- tt„_^ 

^ (~5piO-l) • • • (~5/3n-2Cr„_2) 

.* ^" Hgfi,fe)---V^""-ngt„„2,gt„-J 



X [ku at ■ 



+ kd a; 



d "6 • 



tbl--- tbr,-l 

V-^ iqt^qb,)...^-' 

tbl--- ib„-i 



(36) 



(-1)" kl\at\ 



, ,2, i2^(9bi'9f'2) •••^(%„_2.%„_i; 

+ k^\ab\ 



; y{Qti , qt^ )^ ---y{Qtn-2,it„-i? 

^2 4.2 

tl ' ' ' tn-l 

, y{qti ,qt2)- viQbi ,qb2)--- yiqtn-2 > qtn-i) ■ ^(%. 



f2 ^2 

bl ' ' ' bn-1 



+2kukd{atal + a^ab) 



'n-2 1 qbn-1 



ttl - - - tt„_i tbj_ - - - tb 

n-1 



where ta and tbi are the analogues of equation ([3]) for top and bottom line emission respectively. 
What is particularly appealing about this form is that there are only ever 3 terms, independent 
of the number of final jets. The results are shown in Fig. [12] for Z + 3j and Z + 4j. 



3.4 Higgs Boson Production in Association with Jets 

In this study we will consider just the production of a Higgs boson in-between (in rapidity) jets, 
since this will nicely illustrate the increased predictive power of the factorised current formalism 
over the formalism based on kinematic limits of amplitudes. The production of a Higgs boson 
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(a) (b) 

Figure 10: A comparison of the matrix elements obtained from Madgraph and from our 
formalism for (a) W + 3j and (b) W + 4j. The jets have rapidities A,0, — A in (a), and 
A, A/3, — A/3, — A in (b). The full momentum configurations are in Appendix [Bj 



Pi 




Figure 11: The two diagrams which are added at amplitude level for a Case 2 process with an 
n jet final state. 
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Figure 12: A comparison of the matrix elements obtained from Madgraph and from our formal- 
ism for (a) Z + 3j and (b) Z + 4j (Case 2). The rapidities are as in Fig. [TOj The full momentum 
configurations are in Appendix [Bl 



"outside" jets could be fully included by obtaining results for "Higgs+gluon" -impact factors (see 
Ref. |15j for details) also within the current formalism. 

Let us start by considering the scattering q~Q~ — )• q^Q^H with the scattering matrix 
element 



1l 



9 



4 



{2\u\b), 



(37) 



where the effective coupling of two off-shell gluons to the Higgs boson field through a top-quark 
triangle is described by the vertex obtained in the infinite top-mass limit (please note however 
that the picture of factorised amplitudes does not rely on taking the infinite top-mass limit) 



V^'''%qi,q2) 



S"^"^^ qi.q2-q?q?. 



(38) 



Here, qi = Pa—Pi, q2 = qi —PH- Again, we will use the sum of all allowed helicity configurations 
as the model for all partonic channels. Let us define the spinor/Lorentz string 



•^qQ-^qHQ 



{quq2) = (1 h\t,\a K) g^'^' V^,.,{qi^q2) <7"^'^ (2 h2\v\b h). 



Setting again 



s: 



hahi,^hih2 
qQ^qHQ 



\\SqQ^qHQ{qi,q2)f = ^ 

ha,ha,hi,,h2 

the helicity and colour summed and averaged square of the matrix element is given by 



(39) 
(40) 



M 



qQ^qHQ 



4 (Ag - 1) 



\\SqQ^qHQ{qi,q2) 



1 

h 



6 TT V 



2 1 



(41) 
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Compared to the 2 — )• 2 current scattering of Eq. (jl2p . an extra Lorentz-tensor has been inserted 
in the spinor/Lorentz structure, contracted with the external currents (Eq. (|39|) ). In the square 
of the amphtude, an extra inverse power of the square of the momentum either side of the Higgs 
boson vertex has appeared. This is the same situation as with the emission of a gluon by use of 
the effective vertex Eq. ()24p . Therefore, the approximation for Higgs production in association 
with three jets wih be based on the physical picture of (high energy) current-current scattering 
with rapidity ordered t-channel couplings of the effective gluon emission vertex (Eq. ()24p ) and 
the g*g*H-vertex of Eq. (j38p . in line with the four-jet formalism. Therefore, for the multi-gluon 
emissions, we will take qi,q2 in the argument of the Higgs boson vertex to be the t-channel 
momentum of the propagators either side of the Higgs boson-vertex. 

The approximation for the colour and helicity summed and averaged square of the matrix 
element for the rapidity ordering implied by the subscript 



■^qQ^qHgQ 



\\SqQ^qHQ{qi,q2) 



4 (Ag - 1 



1 / a. \2 1 



ti \6 IT V/ t2 

Ca ^ {-V.V) 1 
1 



(42) 



5' Cf 



where = F^(g2, qs) of Eq. dMI, qi = Pa -Pi, 92 = qi -PH, qs = q2 -P2, and U = qf, with the 
momentum of the quark currents (pa^Pi) and (pb^Pa) respectively. Once again, the factorised 
formalism for the scattering is extremely simple, and the generalisation to multiple emissions is 
straightforward. Again, the approximations for the gluon channels differ only by colour factors 
(multiplication by Ca/Cf for each parton line replaced). 

We show in Fig. [T3]the results of this formalism compared to the full matrix element. We 
also show the approximation of Ref. |18tll9j. We will in Section 14.41 investigate the degree to 
which this formalism captures the physics of the full scattering amplitude. 



3.5 Virtual Corrections and Regular isat ion 

So far, we have discussed how the real emission is approximated to all orders. In this subsection 
we will discuss how the MRK limit of the virtual corrections can be implemented according to 
the Lipatov Ansatz. This will facilitate the construction of a regularised, all-order scattering 
matrix element for each n-parton exclusive final state. We emphasise that:- a) in the current 
study, we construct scattering matrix elements only for rapidity orderings which allow colour 
octet-exchanges between each pair of rapidity ordered particles, b) the contribution from other 
rapidity orderings are systematically suppressed by powers of invariant masses; some of these 
configurations would arise when considering the next-to-leading contributions to the factorisa- 
tion, but can still be included in a t-channel effective description ^17j. The degree to which the 
leading flavour assignments dominate the n-jet cross sections was discussed in e.g. Ref. p^ [T9] 
and are repeated in Section [4.41 We leave the treatment of the sub-leading configurations to a 
later study. 
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Figure 13: The full matrix element for qQ — t- qHgQ, the result of Eq. ()42p ("cur") and the 
approximation of [I8l[l9] ("ADW HE"). 

The Lipatov Ansatz states that the order by order, the virtual corrections to the full n-parton 
scattering amplitude in the MRK limit can be obtained by replacing in the scattering amplitudes 

^ ^ exp[a{qi){yi^i-yi)] (43) 



with 



This ansatz for the exponentiation of the virtual corrections in the appropriate limit of the 
n-parton scattering amplitude has been proved to even the sub-leading level [T6 l [T7 1 [2 H [22] . 

Our discussion of the regularisation of the soft divergences will follow closely the discussion 
in Ref. [19]. We will show that order by order in a^, the soft divergence for the emission 
of gluons cancels with the soft divergence from the virtual corrections implemented according 
to the Lipatov Ansatz for the resummed t-channel propagator. Consider the limit where the 
transverse momentum of the i'th emitted gluon is vanishing. In this limit, 



Pa Pb^Pl ■■■ Pi-1 Pi Pi + 1 ■■■ Pn 



4g^ Ca 



Pa Pb Pi ■■■ Pi-1 Pi + 1 ■■■ Pn 



(46) 



where the matrix element on the RHS has n — 1 final state particles, and p? is the sum of the 
squares of the transverse components of pi in the Euclidean metric. By integrating over the soft 
region p? < of phase space in = 4 + 2e dimensions we find 



1^ 



(2^)2+2.4^ ^'-'•^^'r(l + e)e 

The square of the matrix element on the left hand side of Eq. (|46p contains the exponential 
exp(2a(gj)Ayj_i^j+i). By expanding the exponential to first order in and in e, the resulting 
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pole in e does indeed cancel that of Eq. (j47p . and the combined effect of one soft real emission 
and the first term in the expansion of the Reggeised propagator is a factor 



multiplying the (n — l)-particle matrix element. 

It is clear that the nested rapidity integrals of additional soft radiation in the t-channel fac- 
torised multi-parton amplitudes will build up the exponential needed to cancel the poles from 
the virtual corrections to all orders in Og. The divergence arising from a given real emission is 
therefore cancelled by that arising from the virtual corrections in the Reggeised t-channel prop- 
agator of the matrix element without the relevant real emission. Furthermore, the organisation 
of the cancellation of infra-red poles can be achieved with a simple phase-space slicing. Since the 
t-channel factorised matrix elements are very fast to evaluate and the regularisation procedure 
does not add any complexity (because of the simple IR structure of the t-channel factorised 
matrix elements), the radiative corrections to all orders can be constructed as an explicit phase 
space integral over each number of resolved (|kjp > X'^,i = 2, . . .n — 1) gluons emitted. The 
cancellation of the poles in e ensures that the logarithmic dependence on A generated by the 
lower limit on the transverse momentum integrals cancels with the logarithmic A-dependence 
of the virtual + unresolved-real correction. This is similar to the explicit construction of the 
solution to the BFKL evolution, where the very weak dependence of the solution on A at leading 
logarithmic accuracy was studied in Ref. |23p24]. and in Ref. [25j at next-to- leading logarithmic 
accuracy. 

The construction of an explicit integration over emissions to all orders relies on an efficient 
phase-space generator |181I19]. which should sample final states with the number of particles 
varying by more than an order of magnitude. The problem is significantly different to that of 
a so-called general purpose Monte Carlo (Pythia [1], Herwig [2], SHERPA [3]) where also the 
number of final state particles varies, since in these the approximation to the virtual corrections is 
defined such that the emission of particles is unitary, i.e. does not change the total cross section, 
which allows for a simple probabilistic interpretation. In the problem at hand, (an approximation 
to) the virtual corrections are also calculated, and introduce a suppression of the regularised 
matrix element for any final state with a finite number of partons as the rapidity length of 
the event is increased. This is countered by the (positive) contribution from the emission of 
additional gluons, and introduces a correlation between the number of final state partons and 
the typical rapidity length of an event, which depends also on the transverse momenta. This 
probabilistic correlation is absolutely crucial to incorporate in the phase space generator in order 
to obtain numerical stability in finite amount of time. Such a phase space integrator can be 
efficiently implemented by following the ideas of Ref. [26j . 

4 Applications 

The framework developed here is only useful if it is relevant to the current generation of colliders. 
In this section, we compare the results of using the t-channel factorised matrix elements order by 
order with those obtained with full tree- level matrix elements (extracted from MadGraph |13] ) 
for a 10 TeV proton-proton-collider. We will compare both rapidity and transverse momentum 
distributions for a representative set of processes with a minimal set of cuts. Since the only 
difference between the two calculations presented is the evaluation of the hard scattering matrix 
element, the choice of coupling or parton distribution functions is largely irrelevant for the pur- 
pose of comparison and judging the accuracy of the simple picture for higher order corrections. 




(48) 
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In the later construction of the resummation, we envisage using a NLO set of pdfs. Therefore, 
we choose also in this study to apply the MSTW2008NLO set [27] of pdfs, so that one can judge 
already now the level of corrections which will arise from matching the resummation to fixed 
order accuracy. We will use the fc^-algorithm as implemented in Ref. [25] to define jets. 

In the following, we will only integrate over final states in which the partons are arranged in 
rapidity such that colour octets can be exchanged between all rapidity- neighbours; this constraint 
is relevant only for some partonic channels. Other contributions are systematically suppressed 
by the invariant mass between the non-ordered partons, and the relevance of this suppression 
was checked in an earlier publication [19j. We will repeat this analysis for the updated set of 
cuts in Section [4.41 The non-leading contributions can also be included in a t-channel factorised 
formalism [17J, but we have not yet done so, and so we will compare with the results from full 
tree-level QCD only of the final states which we claim to reproduce. 

We emphasise that the following results will form the basis of a resummation procedure, in 
which the discrepancies between the full tree-level results and those obtained using the approx- 
imation of a t-channel factorised scattering can be corrected for order-by-order, where the full 
result is known. However, the close resemblance between the results obtained order-by-order in 
the simple picture (allowing all-order constructions) and with full QCD should instill trust in 
the results which will be obtained from the resummation procedure. 



4.1 Pure Jets 

Since the two-jet cross section is reproduced exactly for the pure quark channels, we will present 
comparisons only for the three and four-jet channels. We will use the following minimal set of 
cuts 



Pj^ >40GeV \yj\ < 4.5 



In Fig. O we compare the differential cross section with respect to the rapidity difference and 
the azimuthal angle between the most forward and backward jet between the full leading order 
matrix element and the approximation of this obtained in the t-channel factorised approach. 
In the case of ud-scattering (a) — (6), we note that because of the valence quark contribution, 
the rapidity-distribution is peaked at rapidities around 4, and a very good agreement between 
the result obtained within the full tree-level and the simpler t-channel factorised approach; we 
emphasise that no constraints on the minimal separation between jets have been applied. Neither 
a constraint on the similarity of transverse scales in accordance with Eq. ([1]) . 

In the gQ-channel, the t-channel factorised picture obviously works particularly well across all 
rapidities. Also the azimuthal distribution is described very well. Indeed, the (small) discrepancy 
is isolated to the collinear region, where the forward and backward jet is back-to-back in azimuth, 
and thus the third hard jet is radiated in the same direction as one of the other jets. A 
discrepancy in this region is hardly surprising, since the collinear singularity is not included 
in our t-channel factorised description. For the gq and gg-channel, the peak in the rapidity 
distribution obtained with the full LO matrix element is at ever decreasing rapidities; slightly 
above 2 for the channel gg — )• ggg. As expected, in these channels the t-channel dominance is 
only reached for larger rapidities as compared to the ud-channels. However, we emphasise that 
matching corrections can be taken into account in the resummation procedure which will be 
built upon these approximations, such that for smaller rapidities, where the first few orders of 
the perturbative expansion is sufficient (i.e. the effects of the resummation are small), full LO 
results will be used, and for larger rapidities where the resummation effects are important, the 
approximations fare better and so the approximations can be better trusted. This is achieved 
following the matching procedure discussed in Ref. jl9] . 
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Figure 14: A comparison between the differential 3-jet cross section obtained using lowest order 
matrix elements obtained from Madgraph and the t-channel factorised formalism for the process 
(a)-(b) ud — )• ugd with the u-quark incoming in the positive z-direction. (a) the rapidity differ- 
ence between the most forward and backward jet (b) the distribution on the azimuthal angle 
between the most forward and backward jet. Similarly for ug — t- ugg (c)-(d) and gg — t- ggg 
(e)-(f). 
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Figure 15: A comparison between the differential 4-jet cross section obtained using lowest order 
matrix elements obtained from Madgraph and the i-channel factorised formalism for the process 
(a)-(b) ud — )• uggd with the u-quark incoming in the positive z-direction. (a) the rapidity 
difference between the most forward and backward jet (b) the distribution on the azimuthal angle 
between the most forward and backward jet. Similarly for ug — t- uggg (c)-(d) and gg — t- gggg 
(e)-(f). 
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In Fig. [15] we have plotted the similar distributions in the 4-jet cases. We emphasise that 
these (and those in the subsequent sections) are not normalised distributions - the total rates and 
distributions are really approximated by the simple t-channel factorised framework as well as 
indicated in the figures. We note that for all the channels, the peak in the rapidity distributions 
have moved to slightly larger rapidities (roughly one unit) than in the equivalent three-jet 
channels. This is simply because of two effects:- 1) the opening of phase space and 2) that the 
partonic cross section reaches the asymptotic (and large) values at increasing rapidity spans for 
an increasing number of final state particles (see Fig. [2]). We also note that the importance of 
the qg channel compared with the gg-one increases with the jet count. 

4.2 W + Jets 

We now compare the results for the production of a (leptonically decaying) W boson in associa- 
tion with three and four jets. For brevity, we only show plots for the dominant ug — ?■ e'^V(,dg{g)gg 
channel; the t-channel factorised matrix elements reproduce the channels with an initial state 
of only quarks even better. 

We apply the following set of cuts 





> 40 GeV 


Pe± 


> 20 GeV 




< 4.5 


\ye\ 


< 2.5 






Pu± 


> 20 GeV 



The tree-level description of the channels q Q — )• q' Q {W — )•) e are reproduced exactly; 
this is an improvement over the earlier approximations of Ref. [13] , where kinematical limits of 
also the W (and its decay products) had to be applied in order to extract the relevant impact 
factor. In particular, the W had to be assumed to be produced in the forward region, with its 
rapidity increasing with that of the emitting quark. This assumption is in contradiction with 
the requirement of the charged lepton of the W decay products to be central for detection. The 
formalism developed in Section 13.31 requires no kinematical constraints to be placed on the W 
or its decay products in order to extract the building blocks for the resummation. 

We again illustrate the performance of the approximations by analysing the differential cross 
section with respect to the rapidity difference and azimuthal angle between the most forward 
and most backward jets in Fig. [TUJ Also shown are the p± distributions of both the electron and 
the neutrino. The latter is the distribution of missing transverse energy for these events. 

The rapidity distribution is peaked slightly to the right in our formalism, as for the jets, and 
we still see that at larger rapidity differences the results converge. The shapes of the azimuthal 
angle, Pe± and distributions are very good although in each case, our formalism slightly 
underestimates the full matrix element, as was the case for the pure jet events. This is probably 
again due to collinear enhancements in the full matrix element, and the general underestimate 
of the g^-channel for smaller rapidities. 

Figure [T7] repeats these distributions now for the production oi a W boson in association 
with four jets. Again we see that the shapes of the distributions are reproduced quite well. 

4.3 Z + Jets 

We now present similar results for the production of a Z boson (decaying to charged leptons) in 
association with three (and four) jets, again for one of the dominant channels, ug — )• e~^e~ug{g)g, 
with the following set of cuts:- 



Pj± 


> 40 GeV 


Pe± > 20 GeV 


\yj\ 


< 4.5 


\ye\ < 2.5 
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1) 20 40 60 80 100 120 140 160 180 200 20 40 60 80 100 120 140 160 180 200 



P.jtGeV] p^^= Et miss [GeV] 

(c) (d) 

Figure 16: A comparison between the differential W + 3-jet cross section obtained using lowest 
order matrix elements obtained from Madgraph and the t-channel factorised formalism for the 
process ug — ?■ e^vdgg with the u-quark incoming in the positive z-direction. (a) the rapidity 
difference between the most forward and backward jet (b) the distribution on the azimuthal 
angle between the most forward and backward jet (c) the pj_ distribution of the electron and 
(d) the p± distribution of the neutrino, which is the missing transverse energy. 
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kt-jet(R=D.7), pp l^=10TeV 
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Figure 17: A comparison between the differential W + 4-jet cross section obtained using lowest 
order matrix elements obtained from Madgraph and the t-channel factorised formalism for the 
process ug — ?■ e'^vdggg with the u-quark incoming in the positive z-direction. (a) the rapidity 
difference between the most forward and backward jet (b) the distribution on the azimuthal 
angle between the most forward and backward jet (c) the p_L distribution of the electron and 
(d) the p±_ distribution of the neutrino. 
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The differential cross sections with respect to the rapidity difference, the azimuthal angle between 
the most forward and most backward jets and the transverse momentum of both the electron 
and anti-electron are all shown in Fig. [THJ 

The results resemble those for W boson production. The slight difference between the e+ 
and transverse momentum distributions, Fig. 1161 are not seen in Fig. 1181 because identical cuts 
are now applied to both leptons here (the rapidity cut was not applied to the invisible neutrino 
for the W events) . The spin correlation of the W to is also washed out when it is replaced 
by a Z boson. 

Figure [19] repeats these distributions now for the production of a Z boson in association 
with four jets. Again we see that the shapes of the distributions are reproduced very well by 
the simple approximations developed in Section 13.31 

4.4 Higgs Boson + Jets 

In this section, we explore how well our new formalism (section 13. 4p reproduces the results 
obtained using the full matrix element (at tree-level, where the fixed order results can be readily 
obtained). In the following analysis, we apply phase space cuts as follows:- 





> 40 GeV 


Vja ■ Vjb 


< 




< 4.5 


IVja - Vjbl 


> 2 


Vja 


<yh< Vjb 







The merits of these cuts are discussed elsewhere [18l|T3| . In Figure [20] we compare the results 
obtained for the /ij'j-channel in four approximation:- 

1. Lowest order QCD 

2. Lowest order QCD, but including only the flavour and rapidity configurations which are 
taking into account in the t-channel factorised framework 

3. The results obtained using the tree-level (lowest order predictions) of the framework of 
Ref. [I81[I9] 

4. The results obtained using the tree-level (lowest order predictions) of the formalism of 
Section 13.41 

We start by noting that distributions and normalisation is by far dominated by the rapidity 
and flavour configurations of the final states which allow for a colour octet exchange between 
every rapidity-neighbours (these are the only contributions in the leading logarithmic approxi- 
mation). The results obtained by summing over all final state configurations are indicated by 
the full lines, and the contribution from the "leading logarithmic" final states are indicated by 
the dash-dotted line. 

The difference between the results obtained with full QCD and in the t-channel factorised is 
by far dominated by the channels consisting of all gluons plus a Higgs boson. This is consistent 
with the observations made in Figs. [TlVllSI Also, we note that using full tree-level QCD, the 
distribution of the rapidity difference between the most forward and backward hard jet is peaked 
around 2.5 for Higgs boson production in association with two jets, and 3.5 for Higgs boson 
production in association with three jets. The position of the peaks are similar to the situation of 
pure jets, and highlights the universality of the QCD radiation pattern. The rapidity distribution 
obtained with the t-channel factorised matrix element peaks at slightly larger rapidities — this 
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(a) 





Figure 18: A comparison between the differential Z + 3-jet cross section obtained using lowest 
order matrix elements obtained from Madgraph and the t-channel factorised formalism for the 
process ug — ?■ e~^e~ugg with the u-quark incoming in the positive z-direction. (a) the rapidity 
difference between the most forward and backward jet (b) the distribution on the azimuthal 
angle between the most forward and backward jet (c) the p± distribution of the and (d) the 
p± distribution of the e~ . 
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Figure 20: A comparison between the different H+2 jet cross sections using lowest order matrix 
elements obtained from MadGraph, the factorised approximation of [181ll9j and the formalism of 
section [33] for the sum over all relevant subprocesses pp — >• jHj. For the results obtained using 
the full tree-level matrix elements, we have with the line marked "LL conf" indicated the impact 
of taking into account only the flavour and phase space configurations (rapidity orderings) which 
are approximated in the (leading logarithmic) t-channel factorised picture. 
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Figure 21: As in Fig. [20l but now for all subprocesses contributing to pp — )• Hjjj. 
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is simply because the gluon-gluon induced channel is underestimated at smaller rapidities; a 
feature also observed for pure jets. 

We note that the formalism developed in the present study is a slight improvement over 
the earlier study in Ref. |18yi9j in terms of reproducing the shape and normalisation of results 
obtained with full fixed order QCD. It should be noted that the results of the i-channel factorised 
approach are obtained more than 2 orders of magnitude faster than the results relying on the 
evaluation of the full matrix element with MadGraph (all other components of the calculation 
are identical: phase space generation, pdfs etc.). 

5 Conclusions 

We have demonstrated the universal analytical structure of scattering amplitudes in the multi- 
Regge kinematic (MRK) limit of infinite rapidity separation between all produced particles for 
a range of scattering processes which are important for the LHC phenomenology (jets, VF+jets, 
Z+jeis and i?+jets). The universal behaviour in this limit is interpreted in terms of the scat- 
tering amplitudes being dominated by the poles in the t-channel momenta. 

In Section [3] we developed a formalism which is exact in the MRK limit, and fulfils the 
three requirement listed in Section [3] necessary for constructing a relevant all-order summation:- 
Inclusiveness, Simplicity and Accuracy. The formalism is inclusive in the sense that it captures 
the leading real and virtual corrections in the MRK limit. It is sufficiently simple that the 
cancellation of the infra-red poles can be organised explicitly to all orders (because of the simple 
structure of radiative corrections in the MRK limit). The results for the approximations to 
tree-level are obtained more than 2 orders of magnitudes faster than using standard tools for 
evaluating the full tree-level matrix element. This will allow for the all-order sum to be computed 
by explicit evaluation of the exclusive n-parton final states, and thus allow exclusive (or inclusive) 
analyses to be performed. Such an implementation will allow for matching to full tree-level 
results, where such are available. 

We have shown examples of the results obtained in the framework for all 4 processes men- 
tioned above, and compared to the results of using full leading order matrix elements. The 
results are very encouraging; the overall accuracy is good even without phase space cuts to 
enhance the accuracy of the approximation, and the discrepancies which do exists are limited 
to the region of small rapidity separation, where the first few orders in the perturbative series 
should be sufficient and the corrections can be included by matching. This is obviously nec- 
essary in order for the resummation, which will be built on the formalism, to be relevant for 
LHC phenomenology. The phenomenological implications of the resummation will be the focus 
of future studies. 
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A Spinor Representation 



We use the following representation for the spinors. For outgoing particles with 4-momentum 

— E ± pz and p± = Px + ipy, we use 



p, p- 



u~^{p) 



- P. 
p 



'p- 



\ 








and 



u (p) 



I 



( 



\ 








ip- 



lp±l 



For incoming particles with 4-momentum p moving in the + direction, we use: 






and 



u (p) 



\ / 



/ \ 






For incoming particles with 4-momentum p moving in the - direction, we use: 



u^{p) 



'p- 







and 



u (p) 



\ / 

We use the following representation for the gamma matrices: 
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We also use the shorthands 



(ij) = u {pi)u^{pj) and 



u^{pi)u {pj). 
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B Momentum Configurations 

Here we list, for completeness, the momentum configurations used for the plots of the matrix 
elements. Those for the pure jet events (Figs. [5] and [5]) are given in the text. For the production 
of either a W ov Z (decaying to ££) in association with jets (Figs. [3l HI [10] and [T2]). we use the 
following momentum configurations. For 2 jets in the final state : 

Pi = {k±i cosh(yi), k±i cos((/)i), k±i sin(0i), k_i_i sinh(yj)). 



k±i 



40GeV, k^e 



m 



V 



2kj_^{cosh{yg - yf ) - cos{4>^ - (j)e)) ' 
(pi =TT, (t>i = Tr + 0.2, (f>(, = -{tt + 0.2), 
yi = A, 7/2 = -A, = A, ye = A - 1.5, 
P2± = -Pi± - Pi± - Pi±- 



(54) 
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For 3 jets in the final state: 

Pi = ik±i cosh(yj), k±i cos{(j)i), k±i sin((/)i), /c_Li sinh(yi)), 

2 

k±i = k^2 = k^i = 40GeV, k^e, 



2kj_g{cosh{yi - ye) - cos((/)£ - (pe)) ' 
= 27r/3, <^2 = 0, = 7r/2, (j)e = -vr/2, 
yi = A, y2 = 0, = -A, = A, = A - 1.5, 
P3± = -Pi± - P2± - Pi± - Pe±, 
while the 4j final state events use: 

Pi = ik±i cosh(?/i), k±i cos{(j)i),k±i sin((/)j), k±i sinh(yi)), 

m 

2/c_L^-(cosh(y^- - yi) - cos(</>^- - ' 
01 = vr, <?!)2 = vr/2, (p^ = -7r/3, (/)^- = 7r/4, 0^ = -7r/4, 
2/1 = A, 1/2 = A/3, ys = -A/3, y^ = -A, y^- = A, y^ = A 
P4± = -Pi± - P2± - P3± - Pii_ - P£±- 



k±i = k±2 = ^±3 = k^i = 40GeV, k±£ = — — — Y- 



(56) 



For the production of a Higgs boson in association with jets (Figs.[5land ll3p we use the following 
momentum configurations. For 2 jets in the final state: 

pi = (40\/2cosh(A),-40,40,40\/2sinh(A), 

PH = (^402 + m2^, 0, -40, 0) GeV, (57) 
P2 = (40cosh(-A),40,0,40sinh(-A)) GeV, 
and for 3 jets in the final state: 

pi = (40cosh(A),-40,0,40sinh(A), 



PH = (y 402 + mj^ cosh(A/3), 0, -40, ^402 + mj^ sinh(A/3)) GeV, ^^^^ 

P2 = (40cosh(-A/3),0,40,40sinh(-A/3)) GeV, 
P3 = (40cosh(-A),40,0,40sinh(-A)) GeV. 
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